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Abstract 

We consider the formal prolate spheroid differential operator on 
a finite symmetric interval and describe all its self-adjoint boundary 
conditions. Only one of these boundary conditions corresponds to a 
self-operator differential operator which commutes with the Fourier 
operator truncated on the considered finite symmetric interval. 

4 Self-adjoint boundary conditions for 
the prolate spheroid differential operator. 

The study of the spectral theory of the Fourier operator restricted on 
a finite symmetric interval [—a, a]: 

a 

{•S^Exm = ^J e^'^xiO di, teE, E= [-a, a], 



3^e: L\E)^L\E), (4.1) 

is closely related to study of the differential operator generated by the 
differential expression (or formal differential operator) L: 
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The operator L is said to be the prolate spheroid differential operator. 

The relationship between the spectral theory theories of the integral 
operator S^^S^e, E = [—a, a], and the prolate spheroid differential op- 
erator was discovered in the series of remarkable papers jSlPol . [LaPl| . 
[LaP2j . where this relationship has been ingeniously used for develop- 
ing the spectral theory of the operator 3^*^yE- (See also |S12| . |S13| .) 
Actually the reasoning of [SlPo| . [LaPl| . [LaP2j can be easily applied 
to the spectral theory of the operator itself rather the operator 



It should be emphasized that what was used in |SlPo| . [LaPl| . 
[LaP2j this is a certain system of eigenfunctions related to the dif- 
ferential expression L, (|4.2p . These eigenfunctions are known as pro- 
late wave functions. The prolate wave functions themselves were used 
much before the series of the papers [SlPoj . [LaPlj . [LaP2j was pub- 
lished. These functions naturally appears by separation of variables in 
the Laplace equation in spheroidal coordinates. However this was the 
work [SlPoj ■ [LaPlj . [LaP2| where the prolate functions were first used 
for solving the spectral problem related to the Fourier analysis on a 
finite symmetric interval. Until now, there is no clear understanding 
why the approach used in [SlPo| . [LaPlj . [LaP2j works. This is a a 
lucky accident which still waits for its explanation. (See [S13j .) 

Actually eigenfunctions are related not to the the differential ex- 
pression itself hut to a certain differential operator generated by the 
differential expression. This differential operator is generated not only 
by the differential expression but also by certain boundary condition. 
In the case E = (—00,00), the differential operator generated by the 
differential expression — ^ + on the class smooth finite functions 
(or the class of smooth fast decaying functions) is essentially selfad- 
joint: the closure of this operator is a selfadjoint operator. Thus in the 
case E = (—00, 00) there is no need to discuss the boundary condition 
because there is no such boundary conditions. 

In contrast to the case E = (—00, 00), in the case E = [—a, a], < 



o < 00 the minimal differential operator — ^-(l 9 ) ~} 1" is 



symmetric but is not self-adjoint. This minimal operator admits the 
family of self-adjoint extensions. Each of this selfadjoint extensions is 
described by a certain boundary conditions at the end points of the 
interval [—a, a]. The set of all such extensions can be parameterized 
by the set of all 2 x 2 unitary matrices. 

It turns out that only one of these extensions commutes with the 
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truncated Fourier operator S'e, E = [—a, a]. To our best knowledge, 
until now no attention was paid to this aspect. In the present paper, we 
in particular investigate the question which extensions of the minimal 
differential operator generated by L, (|4.2p . commute with L. 



Analysis of solutions of the equation Lx = Xx near singular points. 

For the differential equation 

^ ^_^')^ + t2^(t) =Ax(t), tec, (4.3) 



dt\ a? J dt 

considered in complex plane, the points —a and a are the regular sin- 
gular point. Let us investigate the asymptotic behavior of solutions of 
this equation near these points. (Actually we need to know this be- 
havior only for real t G (—a, a), but it is much easier to investigate this 
question using some knowledge from the analytic theory of differential 
equation.) Concerning the analytic theory of differential equation see 
[Sml Chapter 5]. 

Let us outline an analysis of solution of the equation near the point 
t = —a. Change of variable 

t = —a + s, x{—a + s) = y{s) 

reduces the equation (|4.3p to the form 

d'^v(s) , , dy(s) , , s . ^ 

s-^+P{s)^ + q{s,X)v{s) = ^. (4.4) 

where p{s) and q{s) are functions holomorphic within the disc |s| < 2a, 
moreover p(0) = 1: 

oo oo 
p{s) = l + Y^pus\ q{s) = Y.qk{X)s\ (4.5) 

k=l k=0 

An explicit calculation with power series give: 

1 \a a? A 3 2 

Pi = -^; '?o = y-y, 9i=4 + 4« • (4.6) 

Now we turn to the analytic theory of differential equations. The 
results of this theory which we need are presented for example in [Sml 
Chapter 5], see especially section 98 there. We seek the solution of 
the equation (|4.4p - (|4.5p in the form 

oo 

y{s) = sfJ2 (^kS^ ■ 
k=0 
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Substituting this to the left-hand side of the equation (|4.4p - (|4.5p and 
equating the coefficients of hke powers of s to zero we obtain the equa- 
tions for the determination of p and c^. In particular, the equation 
corresponding to the power is of the form: 

CO = . 

The coefHcient cq plays the role of a normalizing constant, and we may 
take 

CO = 1 . (4.7) 
Equation for p, the so called characteristic equation, is of the form 

= 0. (4.8) 

This equation has the root p = 0, and this root is multiple. According 
to general theory, the equation (|4.4|) - (I4.5I) has two solutions yi(s) and 
2/2(5) possessing the properties: 

The solution yi{s) is a function holomorphic is the disc \s\ < 2a 
satisfying the normalizing condition yi(0) = 1. The solution y2{s) is of 
the form 2/2 (s) = yi{s) In s + z{s), where z{s) is a function holomorphic 
in the disc |s| < 2a and satisfying the condition z{0) = 0. WE may 
calculate explicitly several first coefficients of power expansions 



yi(s) = 1 + ^ Cfcs'', z{s) = ^dk 

k=l k=l 

Aa , - si 

ci = , di = Aa — a H . 

2 2' 2a 

Returning to the variable t = —a + s, we get the following result: 



Lemma 4.1. Let L be the differential expression defined by (|4.2p . and 
X £ C be arbitrary fixed. 

1. There exist two solutions x|f(t,A) and (t, A) of 
the equation Lx{t) = Xx{t) possessing the properties: 

a. The function x^{t, A) is holomorphic in the disc \t + a\ < 2a, 
and satisfy the normalizing condition x^{—a, A) = 1 ; 

b. The function X2 (t, A) is of the form 

X2 (i , A) = Xi {t, A) In {t + a) + w~ (t , A) , 

where the function w^{t, A) is holomorphic in the disc \t+a\ < 2a 
and satisfy the condition {—a, A) = . 
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2. There exist two solutions x^(t,A) and x^(t,A) of 
the equation Lx{t) = Xx{t) possessing the properties: 

a. The function xf (t, X) is holomorphic in the disc \t — a\ < 2a, and 
satisfy the normalizing condition xf[a, A) = 1 ; 

b. The function A) is of the form 

x^{t, A) = xt{t, A) In {a-t)+ w+{t, A), 

where the function w~^{t, A) is holomorphic in the disc \t+a\ < 2a 
and satisfy the condition w~^{a, A) = . 

Given a fixed A, the solutions j;^(t,A), (t,A) are linearly inde- 
pendent, therefore arbitrary solution x{t,X) of the equation (14. 3p can 
be expanded into a linear combination 

x{t, A) = Ci x^ {t, A) + C2X2 {t, A). (4.9a) 

The solutions xf{t,X), x^(t, A) also are linearly independent, and the 
solution x{t, A) can be also expanded into the other linear combination 

x{t, A) = c^x^(t. A) + o^x^it, A). (4.9b) 

Here of, are constants (with respect to t). The solution rr^f (t, A) is 
bounded and the solution X2 (t, A) grows logarithmically as t — > — a. 
Therefore the solution x{t,X) is square integrable near the point t = 
—a. For the same reason, the the solution x{t, A) is square integrable 
near the point t = a. Thus we prove the following result. 

Lemma 4.2. Given A € C, then every solution x{t, A) of the equation 
()4.3|) satisfy the condition 



a 

J \x{t,X)fdt < 00. (4.10) 



Differential operators related to the differential expression L, (|4.11l) . 



With the differential expression (or, in other words, the formal 
differential operator) L, 

d / t^\d 

various differential operators may be related according to whether 
boundary conditions are posed on functions from their domains of 
definition. 



i = -9l(l-^)^ + *'> (4-11) 
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Definition 4.1. The set A is the set of complex-valued functions x{t) 
defined on the open interval {—a, a) and satisfied the following condi- 
tions: 

dx(t) 

1. The derivative of the function x{t) exists at every point t of 

the interval {—a, a); 

dx{t) 

2. The function — - — is absolutely continuous on every compact 

at 

subinterval of the interval {—a, a); 

Definition 4.2. The set A is the set of complex-valued functions x{t) 
defined on the open interval {—a, a) and satisfied the following condi- 
tions: 

1. The function x{t) belongs to the set A defined above; 

2. The support suppx of the function x{t) is a compact subset of the 
open interval {—a, a): (suppx)(s {—a, a). 

Definition 4.3. The differential operator vCmax is defined as follows: 

1. The domain of definition T^n^^^of the operator £max is: 

^Anax = {x ■■ x{t) E L^{{-a,a))r\A and {Lx){t) G L^{{-a,a))}, 

(4.12a) 

where {Lx){t) is defined\jby (14. 2p . 

2. The action of the operator £max is: 

ForxePz:^, £_x = Lx. (4.12b) 

The operator >Cmax is said to be the maximal differential operator gen- 
erated by the differential expression L, (I4.11|l . 

The minimal differential operator C . is the restriction of the max- 

^ mm 

imal differential operator C^^^^^ on the set of functions which is some 
sense vanish at the endpoint of the interval (—a, a). The precise defi- 
nition is presented below. 

Definition 4.4. The operator £^.^ is the closure^of the operator C: 

£„,„ = clos(£), (4.13a) 
where the operator C is the restriction of the operator C^^^ : 

CcCm^,, £ = £i„ax|p^, V^ = V^^^^nA. (4.13b) 



^ Since x e A, the expression {Lx){t) is well defined. 
^ Since the operator C is symmetric, it is closable. 
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By ( , ) we denote the standard scalar product in a, a)): 



For u, V £ L ((—a, a)), {u,v) = / u{t)v{t)dt 



The properties of the operators C^:^^ and /3„jax' 

1. The operator C^^^.^^ is symmetric: 

{C^.^^x,y) = {x,C^.^^y), Vx^y^V^^^j (4.14) 
In other words, the operator C^:^^ is contained in its adjoint: 

C . C (£ . )* ; 

mm — V mm / " 

2. The operators C^^^.^^ and are mutually adjoint: 

(A,.„)*=^n.a., (>C_)*=/:^i„; (4.15) 

In 1930 John von Neumann, |Neu| . has found a criterion for the 
existence of a self-adjoint extension of a symmetric operator and 
has described all such extensions. This criterion is formulated in terms 
of deficiency indices of the symmetric operator. 

Definition 4.5. Let Aq be an operator in a Hilbert space Sj. We 
assume that the domain of definition Vaq is dense in and that the 
operator Aq is symmetric, that isH 

^0 ^ (^o)*. (4.16) 

For complex number X, consider the subspace 

Mx = ^e {{Ao - XI)Va,) , (4.17a) 

or, what is equivalent, 

Mx = {x e^j: {Ao)*x = \x}. (4.18) 

The dimension diuiMx is constant in the upper half-plane Im A > 
and in the lower half-plane ImA < 0; 

diuiAfx = n-\- , ImA>0, (4.19a) 
dimMx = n-, ImA<0. (4.19b) 



3 The relation (|4T6l) means that T^Aq ^ T^A' and Aqx — AqX^/x G T^Aq 
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The numbers and are said to he the deficiency indices of the 
operator Aq, and the subspace A/a is said to be the deficiency subspace 
corresponding to the value A. 

Theorem (von Neumann). 

1. The densely defined symmetric operator admits self adjoint expan- 
sions is and only if its deficiency indices are equal: 

n+ = n_. (4.20) 

2. Assume that a symmetric operator Aq is closed and its deficiency 
indices are equal. Choose a pair of non-real conjugated complex 
numbers, for example X = i, X = —i. The set of all selfadjoint 
extensions of the operator A is in one-to-one correspondence with 
the set of all unitary operators acting from the deficiency subspace 
Mi into the deficiency subspace M-i. In particular, ifn+ = n_ = 0, 
the operator Aq already is selfadjoint. 

We apply the von Neumann Theorem to the situation where the 
operator >Cmin is taken as the operator Aq. Then the equation 

{Aq)*x = Xx 

takes the form 

that is the difi'erential equation 

-^(l- J)^ + iMt) = Ax(t), fe(-a,a), (4.21) 

under the extra condition x{t) € Lp'{—a,a). In particular, the dimen- 
sion of the deficiency space M\ coincides with the dimension of the 
linear space of the set of solutions of the equation (I4.2ip belongings 
to L'^{—a,a). According to Lemma l42l every solution of the equation 
(142111 belongs to L'^{ —a, a). Thus we prove the following 

Lemma 4.3. For the operator >Cmin; the deficiency indices are: 

n+(£„,i„) = 2, n_(£^in) = 2. (4.22) 

Thus, the operator £min is symmetric, but not selfadjoint. The set 
of all its selfadjoint extensions can by parameterized by the set of all 
2x2 unitary operators acting from the two-dimensional Hilbert space 
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Mi onto two-dimensional Hilbert space AA_j, where M±i are defect sub- 
spaces of the operator >Cmin- 

Selfadjoint extensions of operators and self-orthogonal subspaces. 

J. von Neumann, [Neuj . reduced the construction of a selfadjoint ex- 
tension for a symmetric operator Aq to an equivalent problem of con- 
struction of an unitary extension of an appropriate isometric operator 
- the Caley transform of this symmetric operator. This approach was 
also developed by M. Stone, [St], and then used by many others. 

In some situations, it is much more convenient to use the construc- 
tion of extensions based on the so called boundary forms. Especially 
convenient is the usage of this construction for differential operators. 
The first version of the extension theory based on abstract symmetric 
boundary conditions, was developed by J. W. Calkin, [Calj . Subse- 
quently various versions of the extension theory of symmetric opera- 
tors in terms of abstract boundary conditions were developed. The 
dual problem of the descriptions of extensions of symmetric boundary 
relations was also considered. See [RoBj . |Koch| . [Brj . 

Considering the symmetric operator Aq (|4.16|) acting in a Hilbert 
space i^, we introduce the bilinear form form O 

^{x,y)= i^l^^y)-i^^^ly) ^ n-.v^^xv^^^^c. gm) 

The bilinear form is hermitian: 

n{x,y) =n{y,x), yx,yeV^*^, 

and possesses the property 

Q{x,y) = 0, V X € Va* , y e Vaq ■ 

This property allows to consider the form as a form on the factor- 
space £: 

£ = VaJVao. (4.24) 
We use the same notation for the form induced on the factor space £: 

I 

Definition 4.6. The form (I4.23p . is said to he the boundary form. 
The factor space £ is said to he the boundary space. 
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It turns out that 

dim£^ = n_|_ + n_ , (4.25) 

where and n_ are deficiency indices of the operator and that 
the form is not degenerate on £. The non-degeneracy of the form 
means: 

For every non-zero x ^ £, there exists y € £ such that Q{x, y) ^ . 

(4.26) 

Let 5 be a subspace of the factor space £: 

5 C f . (4.27a) 

We identify S and its preimage with respect to the factor-mapping 
T^A* ~^ T^A*/T^Ao (= ^) and use the same notation S for a subspace in 
£ and its preimage in Vaq- 

Vao^SC Va* ■ (4.27b) 

To every S satisfying (|4.27bj) . an extension of the operator is re- 
lated. We denote this extension by Ag: 

Vas=S, AoQAsQAl. 

The operator (As)* , which is the operator adjoint to the the operator 
^5, is related to the subspace S-^: 

{AsT = As^ , (4.28) 

where 5-*"" is the orthogonal complement of the subspace S with re- 
spect to the hermitiam form fl: 

S-^'' = {x££: 17(x,?/) = Vy G 5}. (4.29) 

In particular we prove the following result: 

Lemma 4.4. The extension As of the symmetric operator Aq is a 
selfadjoint operator: As = {As)* , if and only if the subspace S which 
appears in (I4.27bp possesses the property: 

5 = 5^". (4.30) 

Definition 4.7. The subspace S of the space £ is said to be ri-self- 
complementary if it possess the property ()4.30p . 
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It turns out that self- complementary subspaces exist if end only 
if the form f2, (I4.23bp . has equal numbers of positive and negative 
squares. (Which conditions is equivalent to the condition n+ = n_.) 

Thus, the problem of description of all self-adjoint extension of a 
symmetric operator Aq can be reformulate as the problem of descrip- 
tion of subspaces of the space E, (|4.24p . which are self- complementary 
with respect to the (non- degenerated) boundary form 0,, (|4.23bll . 

Selfadjoint extensions of symmetric differential operators. The descrip- 
tion of selfadjoint extensions of a symmetric operator becomes es- 
pecially transparent in the case when this symmetric operator is for- 
mally selfadjoint ordinary differential operator, regular or singular. In 
this case the hermitian form 0, (I4.23y) . can be expressed in term of 
boundary conditions of functions from domain of definitions of the op- 
erator ^Q. This justifies the terminology introduced in Definition 14.61 

We illustrate the situation as applied to the case where the symmet- 
ric operator is the minimal differential operator £min generated by 
the formal prolate spheroid differential operator L. Then the adjoint 
operator Aq is the maximal differential operator -Cmax (See Definitions 
14.41 and [4.31 ) The problem of description of selfadjoint differential op- 
erators generated by a given formal differential operator, has the long 
history. See, for example, [Kr| . |Nail Chapter 5]. The book of [DuSch] 
is the storage of wisdom in various aspects of the operator theory, in 
particular is self-adjoint ordinary differential operators. See especially 
Chapter XIII of |DuSchj . 

In principle we may incorporate the question of description of self- 
adjoint boundary condition for the prolate spheroid differential opera- 
tors in one or other of the existing abstract schemes which is devoted 
to such a description in one or other generality. However to adopt our 
question to such a scheme one need to agree the notation, the terminol- 
ogy, etc. This auxiliary work may obscure the presentation. To make 
the presentation more transparent, we prefer to act independently on 
the existing general considerations and to develop what we need from 
the blank page. 

We use the notations 

^2 

pit) = 1 7T, q(t) = t^, —a<t< a. 

In this notation, the formal differential operator L introduced in (14. lip 
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{Lx){t) = (^p(t)M^^ + q{t)x{t), -a<t<a. 
For every x,y G A, 



{Lx{t)) y{t) - x{t) {Ly{t)) = -[x{t),y{t)l -a < t < a , 



where 



[x{t),y{t)] = -p{t)[^y{t)-x{t)^y 
Therefore, for every x,y £ A and for every a, (3 : —a < a < P < a, 

({Lx{t))^-x{t)(jM^))dt=[x,ym-\x,y]{a). (4.31) 



Lemma 4.5. For every x,y £ T^d^^a.^, there exist the limits 

[x,y]-a= hm [x,y](a), [x, y]" = lim [x, y](/3) . (4.32a) 

Proof. Since x{t), y{t), {Lx{t)), {Ly{t)) belongs to L^((— a,a)), then 

a 

/ {{Lx{t))y{t) - x{t) {Ly{t)) dt < oo. Therefore 
[i{Lxit))W) - x{t) (^Jit)))dt = 



= a^^a+oj - (^2/(t)) j dt. 

Concerning this and related result see for example |HuPy[ Chapter 
10]. □ 

The boundary form Q, constructed from the operator = ^min 
according to (I4.23h|) . 

nLix,y) = ^^max^.^)-(^.-^maxj/) 
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can be expressed in the term of the generahzed boundary values: 

According to (|4.25p and Lemma 14.31 the dimension of the boundary 

space £l: £l = ^£n.ax/^£n.in is: 

dimfi = 4. (4.33) 

To make calculation explicit, we choose a special basis in the space £l 
in which the bilinear form $7/, is reduced to "sum of squares". The 
asymptotic behavior of solutions of the equation Lx = near the 
endpoints of the interval (—a, a), described in Lemma I4T1 prompts us 
the choice of such a basis. Let (p-{t),il}-{t)-,^+{t))4'+{t) be smooth 
functions such that 

= 1, -a < t < -a/2, (/?_(t) = 0, a/2<t<a, 

V'-(t) = ln(a + i), -a < t < -a/2, V'-(*) = 0, a/2<t<a, 

(f+{t) = 0, -a<t<-a/2, (p^{t) = l, a/2<t<a, 

tp+{t) = 0, -a<t < -a/2 , ^l^+{t) = ln(a - t), a/2 <t <a. 

(4.34) 

It is cleat that if x is an arbitrary smooth real valued function, 
then 0/,(x,x) = 0)- In particular, 

^-l(X)X)=0, if X is one of the functions v9_ , V'- ,¥'+, "0+ • (4.35a) 
It is clear that 

^l(X-)X+) = 0) if X± is one of the functions ^±-,'4^± ■ (4.35b) 
Direct calculation shows that 

Ol((^_,'0-) = --, f^L(v?+,V+) = --• (4.35c) 
a a 

Thus, the Gram matrix (with respect to the hermitian form of the 
vectors il)-, ifj^, ■0+ is: 



a 
2 



9.L{^-.,ip-) Ql{^-,'4'-) O.L{ip^,ip+) r2L(99_,V'+) 

^L{iJ-,ip-) Qii'ip-^ip-) QL{'>P_,ip+) VLl{iIj„,^Ij+) 

Ol(vj+,99„) J7l(99+,?/^_) Q.L{ip+,^+) S7l(¥?+,?/'+) 

^l('0+,¥'-) ^^l(V'+,0-) ^L{'4'+,'f>+) ^L{i^+,4^+) 



= J, 

(4.36) 
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where 



J 











i 




(4.37) 



The rank of the Gram matrix is is equal to the dimension of the 
space £l: 

rank J = dimSL = 4 . 

Therefore, the vectors ip-, ip-, V+ generate the space £l = 
^^max/^-Cmin- In particular, the domain of definition T>c^-^^ of the 
minimal differential operator >Cmin can be characterized by means of 
the boundary conditions: 



^£n.in = ■■ <t) G and 

[x, = 0, [x, V'-l-a = 0, [x, ^+]'^ = 0, [x, ^+]'' = 0} , (4.38) 

where the forms [ , [ , were introduced in (|4.5p . 

Lemma 4.6. Let be a bilinear form in the space £ defined by (|4.5p . 
and J be the matrix (|4.37p . 

The vector x^ = (3^ip- + + ^ £l is - 

orthogonal to the vector x"^ = a^(/?_ + Plip- + + P+'>P+ £ "fL, 

i/iai is 

17l(x\x^)=0, (4.39a) 

i/ and only if the vector-row v^i = [ai , /5i , a^, G V is J -orthogonal 
to the vector-row v^2 = [a^, /?^, a^, S V, that is 

v^iJv*2 =0, (4.39b) 

where V is the space of vector-rows equipped by the standard her- 
mitian metric, and the star * is the Hermitian conjugation. 

Thus, the problem of description of self- complementary extensions 
of the operator £jnin is equivalent to the problem of description of O^- 
self-complementarylll subspaces in £, which in its turn is equivalent 
to the problem of description of J-self-complementary subspaces in 

"'As soon as the notion of J-orthogonality of two vectors is introduced, (|4.39bp . the 
notions of J-orthogonal complement and J-self-orthogonal subspaces can be introduced 
as well. 
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C^. The last problem is a problem of the indefinite linear algebra and 
admits an explicit solutions. We set 

Pa 



\{I + J), P- = \{I-J) 



(4.40a) 



More explicitly, 



1 i 

-i 1 

1 i 

-i 1 



P_ 



1 -i 
i 1 
1 -i 



(4.40b) 



i 1 
The matrix J, (|4.37|1 . possesses the properties 

J = J*, j2 = I. 

Therefore the matrices P+, P_, (I4.40al) . possess the properties 

Pi = P+, P^ = P_, P+ = PI, P_ = PI, (4.41) 



P+P_=0, P++P-=I. 



(4.42) 



In other words, the matrices P+, P_ are orthogonal projector matrices. 
These matrices project the space V onto subspaces V+ and V-: 



V+ = VP+, V- = VP_ . 
These subspaces are orthogonally complementary: 

v+ e V- = V . 

The vector rows 



[1, i,0,0], et = [0,0,1, i] 



and 



[1,-^,0,0], 



[0,0,1,-i] 



(4.43) 
(4.44) 
(4.45a) 

(4.45b) 



form orthogonalEI bases in V+ and V_ respectively. 

It turns out that J-self-orthogonal subspaces of the space V are in 
one-to-one correspondence with unitary operators acting from V+ onto 
V_. 



' In the standard scalar product on V = 
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Definition 4.8. Let U be an unitary operator acting from V+ onto 
As the vector-row v runs over the whole subspace V+, the vector 
V + vU runs over a subspace of the space V. This subspace is denoted 
by Su: 

<Su = + vU}, where v runs over the whole V+ . (4.46) 
Lemma 4.7. 

1. Let U be an unitary operator acting from V+ onto V_. Then the 
subspace Sjj is J -self-complementary , that is 

Su = S^-' . 

2. Every J -self-complementary subspace S of the space V is of the 
form Su: 

S = Su 

for some unitary operator C/ : V+ ^ V-| 

3. The correspondence between J -self-complementary subspaces and 
unitary operators acting from V+ onto V- is one-to-one; 

(Ui = U2) ^ {Su, =Su,). 

Proof. 1. The mapping v — > v-\-Uv is one-to-one mapping from V+ onto 
Su- Indeed, this mapping is surjective by definition of the subspace 
Su- This mapping is also injective. The equality v -\- Uv = Q implies 
that f = C/f = since[§ f_L Uv. In particular, dmiSu = dim V+ (= 2). 

If vi and V2 are two arbitrary vectors from V+, then the vectors 
wi = vi -\- viU and 102 = V2 + V2U are J-orthogonal: W1JW2 = 0. 
Indeed, since J = P-f — P_ and = VkP+,VkU = vjJJP^ = 1,2, 
then, using the properties (|4.41l) of P+ and we oobtain 

WiJw*2 = {V^P+ + v'^UP-){P+ - P-)(P>2 + P*-U*V*2) = 
= V1V2 — ViUU*V2 - 

Since the unitary operator U preserves the scalar product, then V1V2 = 
viUU*vl, hence wiJwl = 0. Thus, Su C {Su)^^ - (The symbol _Lj 
means J-orthogonal complement.) Since the Hermitian form (^1,^2) — ^ 
V1JV2 is non-degenerate on V, then dim(5^'') = dim V — dim.Su- Be- 
cause dim V — dim 5(7 = dimSu, we have dim5{/ = dim(5^-^). Hence, 
Su = {Su)'^-' ^ i.e. the subspace Su is J-self-complementary. 

•^Recall that v e V+, C/u e V_, and V+_LV_. 
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2. Let 5 be a J-self-orthogonal subspace. If 



V £ S , V = + V- , v± e V± 



then the condition V-Ljv = 0, that is the condition vJv* = means 
that vivl = V2V2. Therefore, if vi = 0, then also v = 0. This means 
that the projection mapping v — > vP+, considered as a mapping from 
5 — > V+, is injective. For J-self-orthogonal subspace S of the space 
V, the equality dim 5 = dim V — dim 5 holds. Hence dim 5 = dim V+. 
Therefore, the injective linear mapping v — > P+ is surjective. The 
inverse mapping is defined on the whole subspace V+ and can by pre- 
sented in the form v = vi + viU , where C/ is a linear operator acting 
from V_)_ into V_ . This mapping vi ^ vi + viU maps the subspace V+ 
onto the subspace S. 

Since vJv* = 0, then vivl = V2V2, where V2 = viU . Since vi G 
V+ is arbitrary, this means that the operator U is isometric. Since 
dimV+ = dirnV-, the operator U is unitary. Thus, the originally 
given J-self-complementary subspace S is of the form Su, where U is 
an unitary operator acting from V+ to V_ . 

The coincidence Sjjx = <Su2 means that every vector of the form 
vi + V1U1, where vi G V+ can also be presented in the form V2 + V2U12 
with some £ V+: 



Since ui, V2 G V+, fiC/i, V1U2 G V_, then vi = V2, and viUi = viU2- 
The equality viUi = V1U2 for every vi G V+ means that Ui = U2. 



Choosing the orthogonal bases (|4.45ll in the subspaces V+ and V+, 
we represent an unitary operator U by the appropriate unitary matrix: 



The following result is a reformulation of Lemma 14.71 

Lemma 4.8. LetV be the space of four vector-rows, J be a matrix 
of the form (I4.37p . With every 2x2 matrix U = \\upq\\i<p^q<2 , we 
associate the pair of vectors v^{U), v^iJJ): 



Vl + ViUl = V2+ V2U2 . 



Thus, {Su,=Su,)^{Ui = U2). 



□ 



e]_U = ei.nii -|- e'^_U2i, 
e\U = e\ui2 + e^^U22- 



v^{U) = e\+ e\ui2 + e^_U22, 



(4.47a) 
(4.47b) 
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where e^, k = 1,2, are the vector-rows of the form (|4.45|) . and the 
subspace Su ofV which is the linear hull of the vectors v^(U), v^{U), 

Su = \m\\{v^{U), v^{U)). 

1. If the matrix U is unitary, then the vectors v^{U), v'^{U) 
are linearly independent, and the subspace Su is J-self- 
complementary. 

2. Lets be a J -self-complementary subspace of the space V. Then 
S = Sir for some an unitary matrix U . 

3. For unitary matrices Ui, U2, 

{Su, =Su,)^ {Ui = U2) . 
The "coordinate" form of the vectors v^iU), v'^{U) is: 

v^{U) =[l + uu, i{l- uii), U21 , -iu2i ], 

v'^{U)=[ ui2 , -mi2 ,1 + ^22,^(1-^22)]. (4.48) 

Remembering, see Lemma [46l how J-self-complementary subspaces of 
the space V are related to Oi-self-complementary subspaces the space 
£l = ^£max/^.Cmin wc formulatc the following result 

Lemma 4.9. Let us associate with every 2x2 matrix U = ||upg||i<p^g<2 
the pair of vectors d^{U), d'^{U) of the space £l: 

d^{U) = (1 + uii)ip^ + i(l - uii)V'- + U2i(p+ - iu2iip+ , (4.49a) 
d^{U) = Ui2(p- - iUi2lp^ + (1 + U22)f+ + ^(1 - '"22^'+ , (4.49b) 

where the functions ip±, ip± are defined in (|4.34l) . The subspace Qu of 
the space £l is defined as the linear hull of the vectors d^{U), d'^{U): 

gu = )iVi\\{d^{U),d^{U)). (4.50) 

1. // the matrix U is unitary, then the subspace S = Qu is VLl- 
s elf- complementary. 

2. Let S be a Q.L-self-complementary subspace of the space £l. 
Then S = Gu for some an unitary matrix U . 

3. For unitary matrices Ui, U2, 

{Gu,=Qu,)^{Ui = U2). 
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It is clear that a subspace S ^ £l is an Vt^- self- complementary 
suhspace if and only if its Q.L-orthogonal complement S^^l is an ^II- 
s elf- complementary subspace. The subspace {Si/)'^^l can be described 
as: 

= {x(z£^., n^(x,d\u)) = 0, nLix,d\u)) = 0} , 

where d^, are defined in (|4.49l) . (|4.34D . Thus Lemma [49l can be 
reformulated in the following way: 

Lemma 4.10. Let us associate the pair of vectors d^{U), d^{U) with 
every 2x2 matrix U = ||tipg||i<p,g<2 by (|4.49l) . (I4.34p . The subspace 
Ojj is defined as 



Ou = {xe£L- nL{x,d\U)) = 0, nL{x,d\U)) = 0} . (4.51) 



1. // the matrix U is unitary, then the subspace S = Ojj is Vl it- 
self- complementary. 

2. Let S be a Q.L-self-complementary subspace of the space 8l- 
Then S = Ou for some an unitary matrix U . 

3. For unitary matrices Ui, U2, 



Thus there is one-to-one correspondence between the set of all 
2x2 unitary matrices U = \\upq\\i<p^g<2 and the set of all f^L-self- 
complementary subspaces S of the space £l = T^Cn^a.^/'^c^in- This 
correspondence is described as 



where Ou is defined in (li^TD . (lOflll . (O^ . 

On the other hand, the subspaces of the space £l = ^Cma.x/'^Cn^in 
which are self- complementary with respect to the Hermitian form 0^, 
(|4.5|) . are in one-to-one correspondence to self-adjoint differential op- 
erators generated by the formal differential operator L, (|4.1ip . Every 
self-adjoint differential operators C generated by the formal differen- 
tial operator L is the restriction of the maximal differential operator 
J~-max, (|4-^P - on the appropriate domain of definition. According to 
Lemma 14.41 as applied to the operators Aq = Cmin, ^0 ~ ^max, the 
domains of definition of a selfadjoint extension C of the operator >Cmin 
are those subspaces S: 



{Ou,=Ou,)^{Ui = U2). 



(4.52) 



*~-min — — 



■max 



(4.53) 
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which are self- complementary with respect to the Hermitian form f^L, 
(|4.5p . According to Lemma [4.101 Qi-self-complementary subspaces S 
can be described by means of the conditions 

S = {x{t) G : nL{x,dHU)) = 0, nL{x,d^{U)) = 0} , (4.54) 

where d^{U), d^{U) are the same that in (|4.49D . (I4.34p . C/ is an unitary 
2x2 matrix. 

The conditions Q.L{x,d^{U)) = 0, Q.L{x,d'^{U)) = may be inter- 
preted as a boundary conditions posed on functions x G T^Cms.,, ■ Let us 
present these conditions in more traditional form. 

Let as introduce the following notations: 

b-a{x) = lim (t + a)^^, ha{x) = lim (t - a)^, 

t-^-a+O dt t— »a-0 dt 

c-a{x)= lim f(t + a)ln(a + t)^^ -a;(t)y (4.55) 

t-^-a+O \ dt J 

Caix)= lim f(t - a)ln(a - t)^^^ - x(t)V 

t^a-O \ dt J 

Remark 4.1. T/je values 6_a(x), c_a(a;) andba{x), Ca{x) may be con- 
sidered as generalized boundary values related to the function x{t) G 
^^max endpoints —a and a of the interval {—a, a). 

Remark 4.2. The solutions {xi, X2 of the equation Lx = Xx, 

which appears in Lemma \4-1\ satisfy the conditions 

b-aiXi) = 0, C-a{xi) = -1; b^a{x2) = 1, C^a{x2) = 0, 

ba{xf)=0, Ca{xf) = -l; baix:^) = l, Ca(x^)=0. 

Lemma 4.11. For x{t) G ^/Zmax; limits (|4.55l) e^sf, are finite, 
and 

%(X IQi 

b-aix) = —nLix,^p~), c^a{x) = —nL{x,ip-), (4.56a) 

baix) = —^Lix,^p+), Ca{x) = —nL{x,ip+), (4.56b) 

where the functions ip±,4'± o.re defined in (|4.34l) . anc? ^Ae /orm Ql is 
defined by ([45l) . 

Proof. The existence of the limits in (I4.56P follows from Lemma 14.51 
applied to the functions x{t) and y{t) = (p±{t) or y{t) = '4'±(t). The 
equalities (|4.56l) can be obtained by the direct computation using the 
explicit expressions (|4.34p for the functions ^p±{t), tp±{t). □ 



20 



Due to (|4.56p . the equality (|4.36|) can be rewritten as 

-1 0' 
10 

-1 ■ 

10 

(4.57) 

Remark 4.3. The characterization (I4.38P of the domain of definition 
T^Cmin minimal operator £-min can be presented as: 

^£n.in = Ht) ■■ x(t) e and 

b^a{x) = 0, C^a{x) = 0, K{x) = 0, Ca{x) = 0} , (4.58) 

In view of (l436]l . the equahties nL{x,d^{U) = 0, nL{x,d'^{U) = 0) 
take the form 

(1 + Uu) b^a{x) - i{l - Uii) C^a{x) + Ui2 ha{x) + iui2 Ca{x) = , 

(4.59a) 

U21 b-a{x) + iU21 C-a{x) + (1 + U22) ba{x) - i{l - U22) Ca{x) = 

(4.59b) 

Remark 4.4. Since the form Qi{x,y) is anti/inear with respect to 
the argument y: ^L{x,^y) = 7lJ7/,(x,y) for fx €z C, the numbers 
which occurs in (I4.49P must be replaced with the numbers —i, i in ap- 
propriate positions in the equality (|4.59p . For the same reason, the 
numbers Upq which occurs in (|4.49l) must be replaced with the numbers 
(I4.59P . However to simplify the notation, we replace the number 
Upq with the number Uqp rather with the numbers u^. This corresponds 
to that in (|4.49|) we use the matrix U* rather than U as a matrix which 
parameterizes the set of all self- orthogonal subspaces. The matrix 
U* is an arbitrary unitary matrix as well the matrix U. 

Definition 4.9. Let U be a 2 x 2 matrix. The operator Cu is defined 
in the following way: 

1. The domain of definition T^c^ of the operator Cu is the set of all 
xit) e ^?£n,a>c which satisfy the conditions (|4.59ap - (|4.59bl) . (|4.55l) . 

2. For X G T^c^, the action of the operator Cjj is: Cjjx = L^^^X- 

Remark 4.5. In view of (I4.59P and (I4.59|) . for any matrix U , 

T^C ■ ^ T^Cu ■ 



6_a(99_) C^a{^-) b-a{^-) Ca{ip-) 

b-aii)-) C-a{lp~) b-a{lp-) Ca(V'-) 

b-a{^+) C-a{(p+) b-.a{^+) Ca{^+) 

&-a(^+) C_a(V'+) 6-a(^+) Ca{^+) 
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Thus for any matrix U , the operator Cu is an extension of the operator 
r ■ ■ 

'-'mm- 

C-mm Q Q -^max ■ (4.60) 

The equalities (I4.59al) which determine the domain of definition of the 
extension Cu can be considered as boundary conditions posed on func- 
tions X £ ('^66 Remark l4.ll ) 

The following Theorem is a reformulation of Lemma 14.101 in the 
language of extensions of operators. 

Theorem 4.1. 

1. If U is an unitary matrix, then the operator Cu is a self adjoint 
differential operator which is an extension of the minimal differ- 
ential operator £min- J^min C Cu C £max • 

2. Every differential operator C which is a self adjoint extension of 
the minimal differential operator Cmm, is of the form C = Cu 
for some unitary matrix U . 

3. For unitary matrices Ui, U2, 

(Ui = U2) <^ {Cu, = Cu,) . 

Commutator of the operator 3^[-a,a] 3nd Cu- 

Let us calculate the difference !I^_ag^\^C-^a_^x — £iiiax9"[-a,a]^ for ^ ^ 
^£ma>.- Notice that ^max^:: G L'^{[-a,a]), so 3'[_a^a]{J~-ma.xx) is defined. 
Since x E L^([— a,o], the function 3^[_a,a]^(*) is smooth on the closed 
interval [—a, a]. (In fact this function is analytic in the whole real axis.) 
All the more, 3^[-a,a]X G ^.Cmax- Thus for x £ ^£nmx) the difference 

3'[-a,a]^ma.xX - >Cmax3"[-a,a] 2; is well defined. 

Assuming that x G 'C'Cn^a.^ that —a < a < jS < a, we integrate 
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by parts twice[l| 
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rv ^ / 



dx(i) 

For X G T^Cme,^, both limits lim (1 — — exist, are finite, 



and 



lim 

lim(l-tVa')^ 



) 

^ &a(a;) . 



(4.62a) 
(4.62b) 



where b-a{x), ba{x) are defined in (|4.55|) and also appear in the bound- 
ary conditions (I4.59p . Since the limits in (|4.62l) are finite, we conclude 
that \x{t)\ = 0(ln(a^ - t^)) as t ^ ±a, |t| < a. All the more, for 

/ +2 X 

(4.63) 



lim 1-^ x(t) = 0. 

Passing to the limit in (|4.61ll and taking into account (|4.63|) and (|4.62p , 
we obtain 



+ (x)e^"* + 6_(x)e~*"* 
^^\^*«) di. (4.64) 



Transforming the last summand of the right hand side of (I4.64p . we 



Like it is done in (2.31) of the manuscript |KaMa) . 
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obtain 



V a? 



-a 

(„„,e^KV..) = ^(-|,e-.^) = 4(«e..)) 

a 2 



—a —a 

= / -«>'-'^ 4 ( (> - 5) 5 / ^K)'-'^ "'^l - / ^-'^ * ■ 

—a ^ —a ' —a 

(4.65) 

Unifying (|4.64l) and (I4.65p . we obtain the equality 

—a ^ ' 

^ ' —a 

(4.66) 

We summarize the above calculation as 

Lemma 4.12. Let 3"[-a,a] Fourier operator truncated on the 

finite symmetric interval [—a, a] . Let >Cmax the maximal differen- 
tial operator with domain of definition "D/Zn^^x generated by the formal 



differential operator L = — —fl — + . (See Definition 14.31 

dt\ ) at 

If X €z T^c^^^, then 3"[_a,a]^ ^ ^.Cmax; ^'^^ equality holds 

(%a,a]^maxa:)(t) - (^max^a.ap) (t) = '^[b+{x)e'''* + 6_ (x) 



•max 

l2 



e 

(4.67) 



24 



Every selfadjoint differential operator generated by the formal dif- 
ferential operator L is a restriction of the maximal differential operator 
/^max on the appropriate domain of definition. According to Theo- 
rem 14.11 the set of such self-adjoint operators coincides with the set 
of operators Cu, where U is an arbitrary 2x2 unitary matrix. The 
domain of definition of the operator Cu is distinguished from the 
domain T>c^^^ by the boundary conditions (I4.59P constructed from U . 
The next theorem answers the question which operators Cu commute 
with the truncated Fourier operator S^^.^ ^jj. 

Theorem 4.2. 

1. IfU = I, where I is 2x2 identity matrix, then the differential op- 
eratoT^Cj commutes with the truncate(^ Fourier operator 3^[_a,a]- 

3'[^a,a]^I X = £/3"[_a,a] X Vx G • (4.68) 

2. // U ^ I , then the operator Cu do not commute with the operator 

3"[-a,a] •■ 

a. There exist vectors x € T>Cu such that 3"[-a,a] ^ ^£(7; so both 
operators 3"[_a,a]^{/ (^^d Cu'3^[-a^a] 0.'''^ applicable to x, but 

3"[-a,a]'C{/X / Cu3'[-a,a]X ] (4.69) 

b. There exist vectors x £ T)cu such that '3^[-a,a]X T^Cu^ so the 
operator Cu'3^[-a,a] even can not be applied to such x. 

Proof. 

1. For U = I, the boundary conditions (I4.59P take the form 

6_a(x) = 0, baix)=0. (4.70) 

Thus, the domain of definition Vcj of the operator £/ is: 

Vc, = {x: x£ Vc^^^, h.a{x) = 0, ba{x) = 0} . (4.71) 

Every smooth function x{t) on (—a, a) which derivative is bounded: 
sup^g(_„ \x'{t)\ < oo, belongs to T>c^^^- Moreover, according to 
()4.55|) . every such a function satisfies the boundary condition (|4.7ip . 
i.e. h^a{x) = Q,ha{x) = 0. Hence every smooth function on {—a, a) 
which derivative is bounded on {—a, a), belongs to domain of definition 

^Ci ^ Cu for U ^ I. 

^ 3'[-q,q] = Fe for E = [-a, a]. 
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Vcj of the operator Ci. In particular, if x S L'^{{—a,a) and y = 
3^[-a,a]X, then y G 'Dcj- Thus for x £ Vcj both summands in the 
expression 3'[_a^a]^ix — Cj9^[_a^a]X are well defined. Since the operator 
Ci is a restriction of the operator Cmax, then 

In view of (14671) and ([4?70]l . the equahty (1468]) holds. 

2. Let U ^ I. Then at least of one value uu — 1 or U22 — 1 differs 
from zero. For definiteness, let un — 1^0. Set 

7= .//"""w x(t) = ^_(t) + 7</'-W+^o(t), (4.72) 

Z(l - till) 

where a;o(t) is a smooth function which support is a compact subset 
of the open interval {—a, a): 

suppxo <s {—a, a) . (4.73) 

The function xq will be chosen later. According to (|4.57|) . (I4.73P and 
the choice of 7, for any choice of Xo{t), the function x{t) from (I4.72P 
satisfy the boundary conditions (I4.59p . Thus, 

x{t) G . (4.74) 

for any choice of xq . Moreover 

6_„(x) = l, baix)=0. (4.75) 

For the function y{t) = (3"(_a a)3;)(t), the boundary conditions (|4.59p 
either hold, or does not hold. This depends on the choice of the func- 
tion Xq. If (|4.59p hold for this y, then a)^^ £ '^Cu the equality 
(|4.67ll can be interpreted as the equality 

(4.76) 

In view of ^M, {^{-a,a)C-UXm - {Cu'J(^-a,a)X){t) / 0. 

Let us show that both of the possibilities 3"(_(j_a)X G Vc^ and 
'3^{-a,a)X T^Cu are realizable. Since the function y{t) = (9"(_a_a)x)(t) 
is smooth on [—a, a], 

h^a{y) = 0, ba{y) = 0, c_a(y) = -y(-a), c_a(y) = -y{a) ■ 
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Thus as applied to the function y, the boundary conditions (|4.59|1 take 
the form 



(l-nii)?/(-a)-'Ui2y(a) = 0, (4.77a) 
U2iy{-a) - {I - U22)y{a) = Q . (4.77b) 

If, using the freedom of choice of the function xo{t) in (|4.72l) . we can 
arbitrary prescribe the values y{—a) and y{a), then we can either sat- 
isfy the boundary conditions (I4.77P (prescribing y{—a) = 0, y{a) = 0), 
or violate them (if uu ^ 1, we prescribe y{—a) = 1, y{a) = 0, if 
uii 7^ 1, we prescribe y{—a) = 0, y{a) = 1.) The reference to Lemma 
below finishes the proof. □ 

Lemma 4.13. Given complex numbers yi andy2, there exists a smooth 
function XQ{t) on [—a, a] which possesses the properties: 

1. suppxo ^ (—«,«) • 

2. yo(-a) = yi, yoia) = y2, where yo = 3'[-a,a]ixo) ■ 

Proof. The evaluations y{—a) and y{—a) are linearly independent lin- 
ear functionals on the space of functions on (—a, a) which are smooth 
and compactly supported: 

a a 

y{-a) = / x{Oe-'< d^, y{a) = 4^ / me'"^ , 



and the functions e~*"^, e*"^ generating these linear functionals are 
linearly independent on any non-empty open subinterval of the interval 
(—a, a). □ 



Properties of the operator £/. As we have established, Theorems 14.11 
and l4.2l the only selfadjoint differential operator which is generated by 
the formal operator L and which commutes with the truncated Fourier 
operator 3"[-a,a] is the operator Cj. Let as discuss properties of the 
operator £/. 

The following lemma gives an alternative definition of the do- 
main T>Ci- 

Lemma 4.14. Let a function x{t) belong to T^CmB.^- Then the function 
x{t) belong to if and only if both limits 

x{—a) = lim x(t), x(a) = lim x{t) . (4.78) 
exist and are finite. 
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Proof. 1. Functions x{t) belonging to possess the properties 

0, (4.79) 



lim (a2-t2).'^^(*) 



t^±a 



dt 



d_ 



d^ = C^ <oo, C > 0. 



(4.80) 



From (I4.80p and the Schwarz inequahty it follows that 

t2 



tl 



dC <C - h , -a<ti<t2<«- 



All the more, 



dxiO 



di 



< c Vt2 - tl . 



We use the last inequality for ti = —a + 0, t2 = t, where —a < t < a. 
Taking into account (|4.79l) . we deduce that 



_ dx{t) 



Analogously, 



dt 



2 ^2^dx{t) 



dt 



< CVa + t, -a<t<0. 



< CVa - t , < t < a. 



From two last inequalities it follows that 



2 .2^dx{t) 



dt 



1 



< —^Va'^ - t^, -a<t<a. 



Finally, from (|4.79|) and (I4.8QP we deduced the inequalities 



dx{t) 



dt 



and 



< 



C 



c 



-a < t < a . 



(4.81) 



t2 



\x{t2) - X{ti)\ < 



di 



-a<ti<t2<a. (4.82) 
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Since / < oo, the function x{t) is uniformly continuous on the 

— a 

interval {—a, a). Therefore the limits (|4.78|) exist end are finite. 

2. According to Lemma l4.11l both limits lim(^_|_(Q_o)(t±a)^^^^ exist 
and are finite. If x ^ T^Ci, them at least one of these limits is not zero. 
If for example Ximt^a-oit — o)^^^^^ / 0, then the function x{t) grows 
logarithmically as t ^ a — 0. □ 

Theorem 4.3. 

1. The self adjoint operator Ci is an operator with discrete spectrum. 

2. The spectrum of the operator Ci is a sequence of positive eigenval- 
ues of multiplicity one which tends to +oo. 

3. The number A is an eigenvalue of the differential operator Cj if 
and only if there exists the non-zero solution e{t, A) of the boundary 
value problem for the differential equation 

with the boundary conditions 

e(— a, A) is finite , e(a,A) is finite. (4.83b) 

This solution e{t, A) is an eigenvector of the operator Cj corre- 
sponding to the eigenvalue A. 

Remark 4.6. The solutions of the boundary value problem (|4.83p are 
known as the prolate spheroidal wave functions. There is a litera- 
ture where these functions are discussed and studied. See for example 
pISt] . [H], [KPS] . [MSch], [SMCLCj . 

If A is a symmetric operator in a Hilbert space which domain of 
definition Va is dense in and M is a bounded selfadjoint operator 
defined everywhere in , then the operators A and B = A + M {Vb = 
Va) are selfadjoint or not simultaneously, and spectra of A and B are 
discrete or not simultaneously. 

We use this fact in the case when = L'^{{—a,a)), A = £/, 
Mx{t) = x{t) — t'^x{t), so the operator i? is a differential operator 
A of the form 
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which domain of definition Pa coincides with the domain of definition 
V-Dc^ of the operator £/. (See (l47T]l and (|4.12ap .l 

Lemma 4.15. Each of the operators Lj and A is non-negative, and 
for every x £ Vcj = V\ the equalities hold: 



{Cix,x) = J 
{Ax, x) = 



dxiO 



1 



dxiO 



dc 



dC + I eixiOl^dC, (4.85) 

—a 
2 ^ 

di + j k(OPde- (4.86) 



Proof. Let —a < a < (3 < a. Integrating by parts we obtain 
/3 



d_ 

d^ 

1 



e\dx{o 



a? J di 

According to the boundary conditions (|4.70ll . 

\ dx{t) 



f{,_e\dx{o dxio 

e=a J V a^J dC ■ d( 



lim ( 1 - ) = , 

t^±(a-o) V J dt 



According to Lemma [4.141 

\x{t)\ = 0(1) as |t| ^ a - 0. 
Passing to the hmit as a — > —a + 0, /3 ^ a — 0, we obtain the equahty 



d_ 
d^ 



e 



dxiO 



dC 



d^, 



for every x E = Va . (4.87) 



Proof of Theorem [Q Let 

03 = {x G Pa : {Ax,Ax)L2(^_a,a) < l} 



□ 



(4.88) 



30 



be a preimage of the unit ball of the space L'^{—a,a) with respect to 
the mapping x Ax. To prove that the spectrum of A is discrete it 
is enough to prove that the set 53 is precompact in a,a). The 
condition (Ax, Ax)/^2(_q^^) < 1 for a function x G Pa means that 



dc < 1 



(4.89) 



In view of (OTD . 



J -a dflV a 



e_\dx{0\ 



Therefore from (|4.89l) it follows that 



di. 



\x{i)?di<l 



(4.90) 



and 



J-a d^\y a^) di ) 



Inequality (l49T]l the inequality (1480]) for C 
the function x satisfy the inequality 



di < 1 



(4.91) 



According to ([482]) . 



t2 



\x{t2) - x{ti)\ < a3/2 



d^ 



-a < ti < t2 < a . 



(4.92) 



Thus, the set of the functions x belonging to 53 is uniformly bounded, 
(|4.90|) . and equicontinuous, (|4.92p . Therefore, the set *B is precompact 
in L^([-a, a]). 

Thus the spectra of the operators A and Cj is discrete, i.e. consists 
of isolated eigenvalues. According to (|4.85|1 . the eigenvalues of the 
operator Cj are positive. If A is an eigenvalue of the operator Cj and 
e{t, A) is an eigenfunction which corresponds to this A, then, since 
e{t, A) G 'Dcj, the function e{t, A) is continuous in t at the points t = a 
and t = —a. (Lemma 14.141 ) 
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Moreover the function e(t, A) is the solution of the differential equa- 
tion Lx = Ax. As any solution of this equation, the function e(t,A) 
is a linear combination of the solutions x^{t,X) and Xg (t, A). (The 
solutions xf{t, X), {t, X) were introduced in Lemma l4?T]l . Prom the 
behavior of the functions e{t, A), x^{t, A), Xg (t, A) by t ^ — o + we 
deduce that the function e(t, A) is proportional to Xi {t, A): 

e{t, A) = C-X^{t, A), C- / is a constant . 

Analogously, 

e{t, A) = C+xf{t, A), C+ 7^ is a constant . 

Thus, up to the proportionality, there is only one eigenfunction corre- 
sponding to the eigenvalue A. □ 

Remark 4.7. Thus if X is an eigenvalue of Cj, then C-x^ (t, X) = 
C-^-xf [t, X) . Since the differential equation Lx = Xx is invariant with 
respect to the change of variable t —t, then the functions e(— t, A), 
e{t, A) lb e(— t, A) are eigenf unctions as well. Since there in only on 
eigenfunction up to proportionality, then either e{t,X) = e(— t, A), or 
e{t, A) = -e{-t, A). Thus, either C+ = C-, or C+ = -C_. 

Remark 4.8. The spectral analysis of the operator A can he done 
explicitly. Its eigenf unctions are essentially the Legandre polynomials, 
the spectrum also can by found explicitely. The property of the spectrum 
of A to be discrete may be derived from this analysis. However we 
prefer to present less explicit but more general reasoning. 
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